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ABSTRACT 

A new statistical method for pinpointing the massive black hole (BH) in the Galactic Center on the IR 
grid is presented and applied to astromctric IR observations of stars close to the BH. This is of interest 
for measuring the IR emission from the BH, in order to constrain accretion models; for solving the orbits 
of stars near the BH, in order to measure the BH mass and to search for general relativistic effects; and 
for detecting the fluctuations of the BH away from the dynamical center of the stellar cluster, in order to 
study the stellar potential. The BH lies on the line connecting the two images of any background source 
it gravitationally lenses, and so the intersection of these lines fixes its position. A combined search for a 
lensing signal and for the BH shows that the most likely point of intersection coincides with the center 
of acceleration of stars orbiting the BH. This statistical detection of lensing by the BH has a random 
probability of ~0.01. It can be verified by deep IR stellar spectroscopy, which will determine whether 
the most likely lensed image pair candidates (listed here) have identical spectra. 

Subject headings: Galaxy: center — gravitational lensing — infrared: stars 



1. INTRODUCTION 

SgrA*, the unusual radio source in the Galactic Center 
(GC), which was long suspected to be the central massive 
black hole (BH) of the Galaxy, is now securely identified 
with the ~ 3 x 10^ M© dark compact mass that was de- 
tected dynamically by IR observations of stars following 
a Keplerian velocity field (Eckart & Genzel 1997; Chez et 
al. 1998). The identification is based on the non-thermal 
spectrum of SgrA* (e.g. Serabyn et al. 1997), on its com- 
pactness (Rogers et al. 1994), on its lack of detected mo- 
tion (Backer & Sramek 1999; Reid et al. 1999) and on 
its position within 0.1" (1" = 0.04 pc in the GC) of the 
dynamical center of the Galaxy (Eckart & Genzel 1997; 
Chez et al. 1998). The nominal position of SgrA* relative 
to the stars was established by aligning the radio and IR 
grids to within 0.03" {la error) using 4 maser giants in the 
inner 15" of the GC (Menten et al. 1997). Recent mea- 
surements of the acceleration vectors of three stars very 
close to the BH provide another dynamical method for lo- 
calizing the BH and constrain the center of acceleration to 
be inside an elongated ~0.03" x 0.06" (^Icr error) region 
offset by 0.05" from the nominal position of SgrA* (Chez 
et al. 2000). Neither of these positions coincide with any 
observed IR point source. The extreme under-luminosity 
of the BH across the spectrum poses a challenge to the- 
oretical accretion models (e.g. Narayan et al. 1998) and 
motivates the continuing effort to detect, localize and mea- 
sure the fiux of the BH in bands other than the radio. 

A precise determination of the BH position is also rele- 
vant for solving the orbits of stars near the BH, which can 
be used to measure the BH mass M, and the distance to 
the GC i?o (Jaroszynski 1999; Sahm & Gould 1999) and 
to search for general relativistic effects (Jaroszynski 1998). 
The BH will not lie exactly at the dynamical center of the 
stellar cluster due to its Brownian motion. A detection 
of this offset, which is expected to be of order few mas 
(Bahcall & Wolf 1976), can be used to estimate M, and 



to probe the stellar potential and the formation history of 
the system. 

Gravitational lensing of background stars by the BH of- 
fers another, independent way for localizing the BH, which 
is free of assumptions about A/, , Rq or the distribution of 
stars around the BH. The BH lies on the line connecting 
the two images of any background source it gravitationally 
lenses, and so the intersection of these lines fixes its posi- 
tion. The most reliable way to identify lensed image pairs 
is by their spectra, which should be identical up to differ- 
ences caused by non-uniform extinction. At present spec- 
tra are available for all stars brighter than K = 11 mag, but 
only for a few stars fainter than K > 13 mag. No spectral 
identification of lensing has been reported yet. However, 
it is possible to use astrometric and photometric data to 
search statistically for the signature of lensing. This is 
motivated by rough estimates that suggest that there are 
several hundreds of distant background stars within 2" of 
SgrA* (Alexander & Loeb 2001). 

A point mass lens is assumed here since the total stellar 
mass around the BH on the scale of interest is very small 
compared to M,. The possibility of one of the images be- 
ing amplified by a star near the BH (Alexander & Loeb 
2001) is neglected. 

2. METHOD 

The two gravitationally lensed images of a background 
source lie on one line with the point mass lens and the 
source, one on either side of the lens. The angular dis- 
tances of the images from the lens, 6*1.2, and that of the 
source, /3, are related by /3 = 0i_2 — 0\/9i^2 {O2 < 
by definition), where 6'e is the Einstein angle, 6*1 = 
4:GM,Dis/DoiDosC^ = -9i92, M, is the mass of the BH 
and Dis, Dqi and Dos are the lens-source, observer-lens 
and observer-source distances, respectively. The images 
are magnified by Ai^2 = |l ~ ^e/^i,2| (^-g- Schneider, 
Ehlers & Falco 1992). It then follows that the ratios of 
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the angular distances, the projected radial and transverse 
velocities and the fluxes of the two images obey the rela- 
tions 



(1) 



where the velocities are measured relative to the lens. 
These constraints can be used to test whether two point 
sources are a pair of lensed images using astrometric and 
photometric data. The challenge lies in detecting a small 
number of true image pairs among the n observed stars 
against the combinatorial background of the (n — l)n/2 
possible pairs in the presence of measurement errors. 

The constraints on the lensed image pair can be ex- 
pressed more robustly in terms of the projections 

where = (±^1,^2) and w stands for the two-component 

1 /2 

vectors similarly constructed from u,-i.2, '!^'a,2 and 2 j 
with the signs chosen so that C^, = 1 for a lensed image 
pair. The definition of utilizes the fact that stellar 
positions are the best measured quantity in the data an- 
alyzed here. In the presence of measurement errors Eq. 2 
is replaced by a score function, 
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where AC^u is the error in C„, due to the measurement 
errors. The smaller the score, the more likely it is that the 
two stars are a lensed image pair. 

It is not known how many lensed image pairs there arc 
among the stars observed near the BH, if any, and so the 
search for the position of the BH and for the lensed im- 
age pairs has to be conducted simultaneously by looking 
for deviations from a random (unlensed) stellar field. This 
proceeds by enumerating over all possible locations for the 
BH around the nominal origin within a specified search 
field, which is very small compared to the entire stellar 
field; recalculating 6, vt and Vr relative to the trial BH lo- 
cation; calculating Sfj for each of the possible pairs («, j); 
sorting the scores; calculating the likelihood L for the best 
n/e scores under the assumption of a random distribution; 
and finally, identifying the position of the BH as that with 
the minim,al likelihood (ML) in random. The choice of 
the best n/e scores minimizes the number of dependent 
pairs in random draws. Dependent pairs, such as 
and {i,k), increase the noise because they cannot both be 
lensed image pairs. Formal confidence limits on the ML 
position are calculated in the same way as for maximal 
likelihood parameter estimation. 

The probability distribution function (PDF) of the best 
n/e scores and their likelihood in unlensed data are es- 
timated directly by Monte-Carlo (MC) simulations. The 
simulations consist of re-shuffling the actual data by ro- 
tating each stellar position and velocity vector by random 
angles and by randomly permuting the fluxes among the 
stars. In addition, the entire data is randomly shifted 
within the search field to avoid biases with respect to the 
nominal origin. This procedure conserves the averaged 
radial density distribution, velocity field and luminosity 
function while randomizing any angular correlations. 



Some of the pairs are filtered out to further suppress 
the combinatorial noise. A fraction of the stars have spec- 
tral type identification and line-of-sight velocity measure- 
ments. Pairs whose stellar types (early or late type) are 
discrepant arc not compared, and so are stars whose line- 
o f-sight velocities Vz ± Ati^ are discrepant by more than 
VA^u^i -|- /S!^Vz2- In addition, only pairs with ^E,min < 
\/ —B\Q2 < ^^E.max a^c considcrcd. This is necessary for 
avoiding pairs whose distance from the center is too small 
relative to the astrometric errors to allow a reliable test 
of co-linearity, or pairs whose Einstein angle clearly ex- 
ceeds the upper limit Qe = 1.75" ± 0.20" for a source 
at infinity, based on estimates of the BH mass. A'/, = 
(3.0 ±0.5) X IO^Mq (Genzel et al. 2000) and the distance 
to the GC, i?o = 8.0 ± 0.5 kpc (Reid 1993). In practice, 
the results do not depend strongly on the exact filtering 
criteria as long as the search field is not too large. 

The random probability for finding a ML smaller than 
that found in the data is Prand = -Pdata(< ML)A^, where 
f'data(< ML) is the probability of any given position in 
the search field to have a likelihood <ML, which is calcu- 
lated from the random PDF, and where N is the size of 
the enumeration over the BH position. The actual value 
of N is much smaller than the search field grid size (here 
100^) because nearby locations are not independent due to 
the smoothing effect of the measurement errors (Fig. 2). 

However, it can be estimated by A^"^ Psim(<ML), the 
MC median of the likelihood probabilities Psim(<ML) in 
randomly re-shuffled fields constructed from the data. In- 
spection of the likelihood plot in Fig. 2 suggests that there 
are of order few x 10 resolution elements in the search field, 
in agreement with the MC estimate A'" = 28. 



3. RESULTS 

The analysis is based on the compilation of Genzel et 
al. (2000), which includes all available astrometric obser- 
vations of stars in the inner 20" of the GC. The compilation 
lists the stellar positions at a common epoch, proper and 
radial velocities, when available (error weighted averages if 
more than one measurement exists), K-band magnitudes 
and spectral types, when known. The positions are the 
best measured quantity. Individual errors are quoted only 
for the velocities. Genzel et al. (2000) estimate a typi- 
cal astrometric error of 10 mas (random and systematic). 
Magnitudes are accurate to better than 0.1 mag for stars 
with K < 13 mag, while fainter magnitudes are listed with- 
out errors to 0.1 mag precision. A magnitude uncertainty 
of AK = 0.1 mag is assumed here for all stars. The sample 
is augmented by 20 stars observed by Ghez et al. (1998) 
that were not included in the compilation because of their 
uncertain velocity measurements, bringing the total to 295 
stars. Of these, 116 stars have a complete set of measure- 
ments (position, proper velocity and magnitude) and are 
used in the analysis. 

The results presented here are based on 10^ MC runs 
with filtering criteria ^E,max = 2.1" (a 2u margin above 
the maximal value of 0e) and 0E,min = 1-0". The latter 
corresponds to a chance co-linearity probability of ~ 0.3 
(hence the importance of using the additional constraints 
on the velocities and fiuxes). Figure 1 compares the ran- 
dom PDF with the approximating polynom that is used 




for calculating the likelihood, with the distribution of the 
best n/e scores in one of the randomly re-shufHed real- 
izations, with the scores at the nominal origin and with 
those at the most likely BH location. Figure 2 shows the 
likelihood across the 0.3" x 0.3" search field. The posi- 
tion of the ML at ARA = -^0.076", ADec = -0.027" 
relative to the nominal origin lies within the error region 
on the center of acceleration (Ghez et al. 2000) and is 
2.7(7 away from the nominal origin. The ML has a ran- 
dom probability of Prand ^0.01. The position of the ML 
is insensitive to moderate changes in 0E,min, ^E.max, or to 
whether the three innermost stars are included in the data 
set or not, or to higher assumed values of AK up to 0.5 
mag, which may better describe the errors in the flux ratio 
given the patchy nature of the extinction in the inner GC 
(Blum, Sellgren & DePoy 1996). However, Prand increases 
when larger errors are assumed. The results are sensitive 
to the values adopted for the velocities. 90 stars of the 
compiled data were observed with the Keck telescope at 
high resolution (> 2 mas) over a baseline of 2 yr (Ghez 
et al. 1998). The ML analysis of the original Keck data 
has a global minimum at the edge of the search field and 
only a weak local minimum near (+0.08", —0.03"). The 
difference is due to the compiled velocities, which are the 
average of the high resolution, short baseline astrometry 
and of lower-resolution, longer baseline astrometry. 

Table 1 lists the ten most-likely image pair candidates. 
At present none of them have spectra for both images. The 
results of simulations with mock stellar fields with proper- 
ties similar to those of the data (Fig. 2) confirm that the 
ML analysis has a high rate of success in finding the po- 
sition of the BH with as few as 3 image pairs in the data. 
The simulated image pairs are usually ranked very close 
to the top of the list, but some highly ranked candidates 
are false identifications. 

4. DISCUSSION AND SUMMARY 

The results of the combined search for a lensing signal 
and for the BH are intriguing because the most likely BH 



position coincides with the center of acceleration and be- 
cause of the low random probability of the lensing signal. 
However, spectroscopic confirmation of the lensed image 
pair candidates is essential since the statistical significance 
of this result is sensitive to uncertainties in the measure- 
ment errors and because the location of the ML depends on 
the values adopted for the stellar velocities. More gener- 
ally, a comprehensive spectroscopic search for sources with 
matching spectra is the best method for detecting lensing 
since it does not assiime any specific mass distribution for 
the lens. Many of the lensed image pair candidates listed in 
Table 1 have small amplifications, A2 = Ai — 1 <^ 1, with 
a 2 mag average difference between the two images, and 
so none of the second images have spectral identifications. 
Until deeper stellar spectroscopy is available, the statis- 
tical analysis could be improved if individual errors were 
made available for all the measured quantities. It should 
be noted that the 4 maser giants used for aligning the ra- 
dio and IR grids are themselves not stationary (the closest 
to the center, IRS7, has a proper motion of 4masyr~^, 
Genzel et al. 2000). Uncorrected, this could lead to small 
systematic shifts between astrometric observations taken 
at different epochs. Such shifts are unimportant for deter- 
mining the dynamical center, but could introduce errors 
in astrometric lensing searches in heterogeneous compila- 
tions, especially for very high resolution data. 

To summarize, astrometric and photometric observa- 
tions of stars very near the BH in the GC can be used 
to pinpoint the BH on the IR grid with as few as 3 lensed 
image pairs in the field to an accuracy comparable with 
that of dynamical methods. An analysis of all available 
astrometric data indicates that the most likely position of 
the BH coincides with the center of acceleration of stars 
orbiting the BH. Lensing by the BH has been detected sta- 
tistically with a random probability of 0.01. This result 
can be verified by deep stellar spectroscopy. 

I am grateful to M. Fukugita, A. Loeb, P. Sackett and 
E. Quataert for helpful discussions. 
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Fig. 2. — Left: A gray scale plot of In ML (shifted to at the maximum) for 116 stars, as function of the shift in the astrometric grid 
(Ax, Aj;) = (ARA. ADec) over the central 0.3" x 0.3" search field. The X in the center is the origin according to the IR/radio alignment 
with its Ifj error circle (Menten et al. 1997). The polygon is the ~1it error region for the center of acceleration (Ghez et al. 2000). The dots 
are the observed IR sources with their 10 mas errors circles (Genzel et al. 2000). The most likely shift at (+0.076", —0.027") is indicated by 
a cross with Icr and 2(t confidence level contours. Right: As in left panel, for a mock field of 110 unlensed stars and 3 lensed image pairs, 
with the BH at (+0.076", -0.027"). 



Table 1 

The ten most likely image pair candidates. 
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"Einstein angle and radial distances in arcseconds, relative to best fit BH position. 
''The amplification of the brighter image. For the second image A2 = Ai — 1. 
'^K-band magnitudes. A spectrum also exists if marked by (f). 
''Radial and transverse velocities in kms~^, relative to best fit BH position. 
''Row number of source in table 1 of Genzel et al. (2000). 
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